Anisotropic Cosmological Models with Energy Density 

Dependent Bulk Viscosity 



Wung-Hong Huang 
Department of Physics 
National Cheng Kung University 
Tainan ,70101, Taiwan 



An analysis is presented of the Bianchi type I cosmological models with a bulk viscosity 
when the universe is filled with the stiff fluid p = e while the viscosity is a power function 
of the energy density, such as rj = a\e\"'. Although the exact solutions are obtainable only 
when the 2n is an integer, the characteristics of evolution can be clarified for the models 
with arbitrary value of n. It is shown that, except for the n = model that has solutions 
with infinite energy density at initial state, the anisotropic solutions that evolve to positive 
Hubble functions in the later stage will begin with Kasner-type curvature singularity and 
zero energy density at finite past for the n > 1 models, and with finite Hubble functions 
and finite negative energy density at infinite past for the n < 1 models. In the course of 
evolution, matters are created and the anisotropics of the universe are smoothed out. At 
the final stage, cosmologies are driven to infinite expansion state, de Sitter space-time, or 
Friedman universe asymptotically. However, the de Sitter space-time is the only attractor 
state for the n < 1/2 models. The solutions that are free of cosmological singularity for any 
finite proper time are singled out. The extension to the higher-dimensional models is also 
discussed. 



E-mail: whhwung@mail.ncku.edu.tw 

Publised in: J. Math. Phys. 31 (1990) 1456-1462 



1 



1 Introduction 



The investigation of relativistic cosmological models usually has the energy momentum ten- 
sor of matter as that due to a perfect fluid. To consider more realistic models one must take 
into account the viscosity mechanisms, which have already attracted the attention of many 
investigators. Misner [1] suggested that strong dissipative due to the neutrino viscosity may 
considerably reduce the anisotropy of the blackbody radiation. Viscosity mechanism in the 
cosmology can explain the anomalously high entropy per baryon in the present universe [2,3]. 
Bulk viscosity associated with the grand-unified-theory phase transition [4] may lead to an 
inflationary scenario [5-7]. 

An exactly soluble isotropic cosmological model of the zero curvature Friedman model 
in the presence of bulk viscosity has been examined by Murphy [8] . The solutions that he 
found exhibit an interesting feature that the big bang type singularity appears in infinite 
past. Exact solutions of the isotropic homogeneous cosmology for the open, closed and flat 
universe have been found by Santos et al .[9], when the bulk viscosity is the power function 
of energy density. However, in some cases, the big bang singularity occurs at finite past. It 
is thus shown that Murphy's conclusion that the introduction of bulk viscosity can avoid the 
initial singularity at finite past is not, in general, valid. (The extensive collections of exact 
isotropic solutions are those found in Ref. 10.) 

Belinskii and Khalatnikov [11] analyzed the Bianchi type I cosmological models under 
the influence of viscosity. They then found the remarkable property that near the initial 
singularity the gravitional field creates matters. Using a certain simplifying assumption, 
Banerjee and Santos [12-13] obtained some exact solutions for the homogeneous anisotropic 
model. Recently Banerjee et al.[14] obtained some Bianchi type I solutions for the case of 
stiff matter by using the assumption that shear viscosity is the power function of the energy 
density. However, the bulk viscosity coefficients adopted in their model are zero or constant. 

In this paper, without introducing the shear viscosity, we shall examine the Bianchi type 
I cosmological models with bulk viscosity {rj) that is a power function of energy density (e), 
i.e., 1] = a|e|", when the universe is filled with the stiff matter p = e. We are interested in 
the cosmological solutions that will eventually go to the states of positive Hubble functions. 
The exact solutions are obtained when n is an integer. Furthermore, through some analyses, 
we can know how evolutions of the models with arbitrary value of n will be. We prove that 
the isotropic de Sitter space-time is a stable attractor state as i ^ oo if n < 1/2. It is thus 
in accord with the "cosmic no hair" theorem [15-17] even though the strong energy condition 
[18] is violated [19]. (The weak energy condition question in the anisotropic viscous models 
has been discussed by Barrow [20].) 

We also find that, for the models oi n < I there are solutions (in fact, all solutions in 
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< n < 1/2 models) that can avoid the cosmological singularity at any finite proper time. 
Our models show that the anisotropies of the universe are smoothed out and matters are 
created by the gravitational field in the course of the evolution, in agreement with the results 
found by others [11,14]. 

The models of n = 0, 1 have been discussed in our previous paper [21]. However, the 
method invented there cannot be used to solve the models with other n, and the analyses 
of the n = model were incomplete. Barrow [22] had also given further discussion of bulk 
viscous models in theories possessing quadratic curvature. It is worth mentioning that the 
n = 3/2 model can be used to describe the quantum production of infinitely thin Witten 
strings [23] on super-horizon scales in the very early universe (see the arguments in the paper 
of Turok [24]), which has been analyzed by Barrow recently [10]. 

The outline of this paper is as follows. We first discuss the axially symmetric Bianchi type 

1 model in which there are ony two cosmic scale functions. In Sec. II we reduce the Einstein's 
field equation to a pair of coupled differential equations that become an integratable equation 
when we define two suitable variables in Sec. III. The exact solutions for the models with 
integer In are then found. The analyses about the initial and final states of the models with 
any n are given in Sec. IV. To discuss the generic Bianchi type I model with multiple cosmic 
scale functions we then in Sec. V present a simple method that can also be used to analyze 
the higher-dimensional models. Section VI is devoted to conclusions. 

2 Einstein's Field Equation and Some Analyses 

We first consider the axially symmetric Bianchi type I model with a metric in the form 



where X and Y are functions of cosmic time t alone. The field equations to be solved are 



dt + X{tfdx^ + Y{tf{dy^ + dz^), 



(2.1) 



((e - p)/2)gfj,u + (e + p)u^UnU, 



(2.2) 



where e is the energy density, and ti^ is the four-velocity that satisfies 



(2.3) 



The total pressure p is defined by 



p = p — rju'. 



(2.4) 
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where p is the pressure coming from the perfect fluid and rj is the bulk viscosity. Choosing 
a comoving frame where ti^ = 5° , the Einstein's field equations (2.2) lead to 

^ + HW^{e-p)/2, (2.5) 

^ + hW = {e-p)/2, (2.6) 

_ (2//2 ^ ^ 2e, (2.7) 



where H and /i are the Hubble functions defined by 

/.^fr- (2.8a) 

h = ^X-\ (2.8&) 

and W is the total expansion function 

W = 2H + h. (2.8c) 

We also have a relation 

e-p^{2--f)e + r)W, (2.9) 



where 7 is defined by the equation of state 

p=(7-l)e, /<7<2. (2.10) 

In this paper we only consider the stiff matter, i.e., 7 = 2, which is the possible relevance 
of the equation of state p = e regarding the matter content of the early universe [25-26]. 
Using the assumption that the bulk viscosity is a power function of energy density, 

77 = ae",n>0. (2.11) 
Equations (2.5) and (2.6) then become 

^ + HW ^^[H{H + 2h)]''W, (2.12) 

^ + hW =^[H(H + 2h)rW. (2.13) 
at 2 



4 



The work that remains is to analyze the above coupled differential equations -and find 
their exact solutions. We will express these solutions as the flows in the phase space H x h 
and thus find the dynamical evolutions of cosmology. We only consider the physical plane 
where the trajectories shall evolve to the positive Hubble functions in the latter stage. The 
regions where e < that violate the weak energy condition[18] are needed for the solutions 
of n < 1 models as discussed in next section. We then find there that the evolutions of the 
cosmology arc confined in the regions oi H > and H + 2h > for the n > 1 models, while 
these for the n < 1 models are > 0. 

Let us first determine the singular points, fixed points, and cosmic time in the phase 
space. 

A. Singular points 

Equation (2.7) tells us that the energy density becomes infinite only if H and/or h are 
infinite. From Eqs. (2.12) and (2.13) we also know that dH/dt and dh/dt can be infinite 
only if H and/or h are infinite. As the Riemann scalar curvature can be written as 

dt 

we sec that R can be infinitely large only if H and/or h arc infinite. Therefore, the singularity 
of diverge R and e could occur only if H and /or h is infinite. 

B. Fixed points 

The fixed points are the solutions of Eqs. (2.12) and (2.13) once we let dH/dt = dh/dt = 
0. They are 

{1)W = => H = h = or 2H=-h^0, (2.14a) 
{2)H ^h^HD = [(372)a]^/(^-^'^), n ^ 1/2, (2.146) 
(3)if — h — any value if a — ac = 2/-\/3, n — 1/2 (2.14c) 

Using Eqs. (2.12) and (2.13) we can determinate the signs of dH/dt and dh/dt in the 
neighborhoods of the de Sitter space-time ol H — h — Hd and thus determine the stabihty 
of the de Sitter state. It is then found that the isotropic cosmologies with n < 1/2 display 
inflationary behavior but those with n > 1/2 shall exhibit the deflationary behavior as found 
in Ref.9. When n — 1/2, then H — h — and H — h — oo is the attractor state for a < ac, 
and a > ac, respectively. The investigation appears to display the division into n > 1/2 , 
n — 1/2, and n < 1/2, also discussed by Barrow [19.20]. 

C. Cosmic time 
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The solutions expressed as the trajectories on the phase plane do not explicitly depend on 
the cosmic time. However, through a simple analysis we can determine whether the proper 
time in a solution is finite or infinite. 

Equation (2.12) can lead to 

rH dff rt 

J [{a/2)[H{H + 2h)Y-H]W^ J "^^^ ^^'^^^ 
We then see that the left-hand side in the above equation will become infinite only if 
is on a fixed point, and one can show that the cosmic time of a point on a trajectory 
corresponding to a solution will be infinite if the trajectory has already met a fixed point. On 
the contrary, the trajectory starting with diverge H will have finite proper time. 

3 Exact Solutions 

We now begin to solve Eqs. (2.12) and (2.13). After dividing the former one by the latter 
one, we obtain 



dH _ a[H{H + 2h)f -2H 
It ~ a[H{H + 2h)Y -2h ■ 

The above equation can lead to 

dL _ ?,a[HLf~2L 
'dt ~ a[HL]'' - 2h ' 
dB _ -2B 
'dt ~ a[HLY-2h' 



(3.1) 



(3.2a) 
(3.26) 



where 

L = H + 2h, B = 2{H-h). (3.3) 
Dividing Eq. (3.2a) by Eq. (3.2b) one gets 



dL_L 
dB~B 



1- -3^-"L2"-^(1 + -)" 
2 L 



(3.4) 



Using a new variable 

A^L/B, (3.5) 

then Eq. (3.4) gives a simple form as 
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dA/{A' + AY = -{a/ 2)3^-'' B'^'-'dB. (3.6) 

The variables are thus separated and the equation becomes integratable. Through the inte- 
gration by part, the exact solutions can be found when 2n is an integer. 

Since we will describe our solutions in the phase plane H x h, we express H and h in the 
variables r and 9: 

H — rcosO, (3.7a) 
h — rsinO, (3.76) 



in terms of which A and B become 

2(1 - tan0) ' ^ ' 

B ^2r{cose - sine). (3.9) 

Therefore, if Eq. (3.6) can be integrated exactly, the solutions that relate function r to 
the variable 9 will be found, and trajectories in phase plane can be plotted exactly, which 
in turn determine the evolutions of the cosmology. Various integration constants chosen in 
integrating Eq. (3.6) will produce the various trajectories which correspond to the solutions 
of the same model but with different initial states. The arrows in the trajectories, which tell 
us the directions of the evolution of the cosmology, are easily determined from Eqs. (2.12) 
and (2.13). 

As argued below, we also need to consider the solutions of negative energy density. A 
natural extension is the models with the viscosity function rj — «(— e)". Through the same 
procedure we can find the following relations: 

dA/{-A^ - AY ^ -{a/2)3'-''B^''-^dB {if e<0;H,W >0), (3.10a) 
dA/{-A^ - AY = (a/2)3^-"S2"-2(^S {if e,H <0;W >0), (3.106) 

The variables arc separated and the exact solutions can be found when n is an integer. 
As examples, we will give some explicit-solutions: 
(1) n = 

h-HD = C{H-HD), (3.11) 
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where C is an integration constant, and is defined in Eq.(2.14b). The solutions are thus 
all the straight lines that pass through the point oi H = h = This means that the initial 
state of cosmology shall begin with H and/or h oo. Therefore, initial singularity will arise. 
However, we must be careful now. The analyses in the above section show that the points 
on the line of = are the fixed points. Therefore some anisotropic solutions will begin 
on these points (in which e is negative) at infinite past; during the evolution the cosmologies 
are isotropized and driven to the de Sitter state asymptotically. Using the relation of Eq. 
(3.11) we can from Eq. (2.12) find the function H{t) which then explicitly shows this fact. 

Note that the energy density of these solutions is negative in the early stage; this will 
violate both the weak energy condition and the strong energy condition, and there is no 
singularity at any finite proper time. 

Although it is difficult for the negative energy density to appear classically, it could be 
found in the quantized matter field. Also, as Hu [27] has discussed, the quantum dissipative 
process of the particle production could be formulated in terms of relativistic imperfect fluid. 
Accordingly, it seem that there will be, more or less, some quantum senses in these solutions. 
It is interesting to mention that the introducing of quantized matter field into the energy 
momentum tensor can sometimes lead to avoidance of the cosmological singularity, as found 
by Parker Fulling [28]. 

(2) n = 1/2 



r = C[±{cos9 - sm^)]("^/")-^[Vcos^ + 2sin9 + sVcosO]'^"''/'^ {if e > 0), (3.12a) 



where C is an integration constant and ac is the constant defined in Eq. (2.14c). The 
solution shows that H and/or h can go to infinity only if H = h, i.e., cos6 = sinO. Therefore, 
this model is with finite Hubble functions at the states of zero energy density {H — Q or 
H + 2h = 0). 

The experiences from the analyses of the n = model tell us that, as the points on 
the lines of zero energy density are not the fixed point, it is now also needed to discuss the 
regions of e < 0. (Note that we only discuss the expanding solutions, i.e., W > 0.) Wc then 
must analyze the extended model in which the bulk viscosity is r] = «(— e)^/^. Equation 
(3.10) can give the explicit solutions: 



C 



r — 



sinO — cosO 



exp 



dc _,,2 + tan9 

± — sin [- -) 

a 1 — tanU 



(3.126) 



where C is an integration constant. 
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After plotting the flows in the phase plane (see Fig. 1, in which, for clarity we adopt the 
non- uniform scale), we then find that anisotropic cosmologies shall always begin with finite 
negative energy density at the states with H + h = in infinite past; the Riemann scalar 
curvature and Hubble functions are finite in the initial phase. During the evolution, the 
energy density is increasing subsequently and the anisotropics of the universe are smoothed 
out. At the final stage as t oo, depending on the value of a (and not on the value of 
C), there arc three classes of states that may be approached asymptotically: 1. Both the 
energy density and Hubble functions go to infinity if a > ac- 2. Energy density is finite 
and space-time is attracted to a de Sitter universe (which is a function of the value of C) 
if q; = ac- 3. Both the energy density and Hubble functions decrease to approach zero and 
the model is driven to the isotropic Friedman universe a < a^. 

The last two cases provide us with the solutions that are free cosmological singularity for 
all finite proper time. 

(3) n = 1 

r = [a{cos9 - sine]-^[C - ln{l + 2tan9)], (3.13) 



where C is an integration constant. After plotting the fiows in the phase plane (see Fig. 
2, where, for clarity, wc adopt the nonuniform scale), we then find that the anisotropic 
cosmologies shall start from the vacuum states and end in another fixed point or infinite 
expansion state. It is then found that, except in the isotropic model {h = H) that has been 
investigated by Murphy [1], the cosmologies shall always begin with zero energy density 
at the initial phase of singularity. During the evolution, the energy density is increasing 
subsequently and the anisotropics of the universe are smoothed out. At the final stage as 
t ^ oo, depending on the integration constant C (and not on the value of a), there are three 
classes of state that may be approached asymptoptically: 



FIG. 1. The phase plane trajectories of the n — 1/2: All solutions will begin with 
W — O,e<oeitt^ — oo. 
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FIG. 2. The phase plane trajectories of the n—1 modeL The solutions that evolve into 
positive Hubble functions at later stage will begin with e = at finite past. 



1. Both the energy density and Hubble functions go to infinity if H > h, C > ln3 or 
h> H,ln3>C>0. 

2. Energy density is finite and space-time is attracted to a de Sitter universe [which is 
determined by Eq. (2.14b) ] if C = ln3. 

3. Both the energy density and Hubble functions decrease to approach zero and model 
is driven to the isotropic Priedmann universe \l H > h, InS >C>Qorh>H,C> In3. 

The solutions that start at fixed point on the line of 1^ = are also described by Eq. 
(13) if one lets a — > —a. However, the trajectories that begin with 1^ = do not go into 
positive Hubble functions states at the final stage. 

(4) n = 3/2 

r = 

a {cos 9 — sinOy 

This model can be used to describe the quantum production of infinitely thin Witten 
strings [23] on super- horizon scales in the very early universe [24,10]. 

(5) n= 2 

r = 

a{cos9 — sinOy 

and so on. 

All the models of n > 1, as will be proved in the next section (for any n), possess the 
same characteristics such as the isotropization of the cosmology, beginning with zero energy 
density and with infinite Riemann scalar curvature, ere ating the matters in the course of 
evolution, and having three classes of final state. 



4(2 + tanO) 
Vl + 2tane 



(3.14) 



3(2 + tane)(l-tane) 9, 

^ ; ^ + - 1 + 2tane) - c 

1 + 2tan9 2 ^ ' 



(3.15) 



4 Analyses of Initial and Final States 

We give in this section the analyses of the initial and final states. The results can clarify the 
characteristics of the cosmological evolutions for the models of any values of n. 

A. Initial states The key equation to analyze the initial and final states of our model is 
Eq. (3.6). It can lead to 

(2A + 1) , ^ {2A + 1) 



^1 + ^yd^ + C2{2A + 1){A^ + AT + ^3^^^^ + C,-^^, ^ 
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^'^'WTW^ + - = CoB'^-' = Co[2r{cos0 - sin9)]'^-\ (4.1) 

where Ci are constant numbers which only depend on n, and a is chosen to satisfy < o" < 1. 
It is important to notice that C3 is nonzero only if n > 1, C4 is nonzero only if n > 2,..., 
and so on. If the universe is in the initial state if + 2/i — > 0, then A — > and + A — > 0, 
and we can easily prove that the first term (neglects the integration constant) and second 
term on the left-hand side of Eq. (4.1) shall approach zeros while the other terms become 
infinite. Therefore the value of r in Eq. (4.1) is finite for the states that have vanish energy 
density, if and only if < n < 1. 

However, it is now necessary to discuss the solutions in the regions of e < 0, as the e = 
state is not the fixed point on the phase plane. Prom Eq. (3.10) one can show that r is finite 
as — > 0. Hence we have proved that the cosmologies of < n < 1 shall begin with finite 
negative energy density and zero total expand function at infinite past . They will then go 
into the positive energy density state. Although the models of n > 1 shall begin with zero 
energy density and diverge values of H and h , one can from Eq. (2.15) conclude that the 
models with n> 1 shall begin at finite past. 

Furthermore, for the ngeq > I models we can from Eqs. (2.5) and (2.6) prove that near 
the initial phase the curvature singularity is the Kasner type, although the energy density is 
zero. This fact was first found in the letter of Belinskii and Khalatnikov [29], in which only 
the anisotropic n = 1 was analyzed. 

B. Final states To analyze the final states we consider three cases separately: 

Case 1: n > 1 The key equation can lead to 

r why + ^^WTW-^ ^ ^3 (If^ ^ - - ■ l^-^) 

where Cj are constant numbers that depend only on n, and a is chosen to satisfy 1 < cr < 2. 
The value of C2 is nonzero only if n > 2, C3 is nonzero only if n > 3,..., and so on. As 
universe approaches isotropic state, i.e., H ^ h, then A ^ 00 and S — > 0, which in turn 
imphes that the first term (neglects the integration constant) on the left-hand side ofEq. 
(4.2) becomes zero. Therefore, depending on the chosen integration constant, r may be 
??oo. This situation is hke that in the model oi n — I. Therefore, depending on the initial 
state, the cosmology may be driven to infinite expansion state, de Sitter space-time, or 
isotropic Friedmann universe at the final stage. 

Case 2: 1/2 < n< 1 

The key equation can lead to 
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dA 



+ C2 



{2A+1] 



(A2 + Ay ^A2 + A) 



= Co[2r{cos9 - sin9)] 



2n-l 



(4.3) 



where Cj, are constant numbers that depend only on n, and a — n + 1, thus 1 < a < 2. 
Using the arguments hke that in Case 1. we can also show that the models with 1 > n > 1/2 
have three classes of final states. Hence we have proved that, depending on the initial state, 
the cosmologies for the models ofn > 1/2 may be driven to infinite expansion state, de Sitter 
spacetime, or isotropic Friedmann universe at the final stage. The model of n = 1/2 was 
discussed in Sec. III. 

Case3: n< 1/2 

The key equation can lead to 



n 



2(n- 1) 



dA 



A (A2 + n(A2 + Ay 



-3 



1— rt 



a 



2(2n- 1) 



B 



2n-l 



(4.4) 



As universe approaches isotropic state, i.e., H — h {it implies A — > 00 and S — > 0), it is easy 
to prove that the first term in the bracket becomes zero (neglects the integration constant) 
and 



[3"2- 



-("+l/2)Q,ll/(l-2n) 



(4.5) 



no matter what the value of integration constant that will be chosen. The state corresponding 
to Eq. (4.5) can be easily checked to be just the fixed point defined in Eq. (2.14b). Hence, we 
have proved that solutions of the models ofn < 1/2 shall always be attracted to an isotropic 
de Sitter state at final stage. 

For the n < 1 models there are the solutions that start with 1^ = at infinite past and 
then are attracted to the original point on the phase plane. However, they never go into the 
state of positive Hubble function. 

5 Models with Multiple Hubble Functions 

The methods described in the above sections can only be used to study the models with 
two Hubble functions. We will now give a simple algorithm that can be used to analyze the 
models with multiple Hubble functions. 

Let us consider the D+l-dimensional Binachi type I models. The Einstein's field equation 
Eq. (2.2) can lead to 
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^ + H,W^{e-p)/2, i^l,2,...D, (2.5) 
W'-j:Hf = e, (5.2) 

i 

where Hi, are the Hubble functions and W is the total expansion function. With the relation 
Eq. (2.9), and letting 7 = 2 in there, we can from Eq. (5.1) find 

dlnjHj - H,) ^ dlnjHj - Hk) 

dt dt ' ^ ' ' 

This equation tells us that one can express all other D — 2 Hubble functions in terms of only 
two Hubble functions. One can see that this is a very general property and may be used to 
analyze many other anisotropical cosmological models As examples, we will now describe the 
procedure of how to determine the three Hubble functions in the four-dimensional model. 
Equation (5.3) can yield a relation 

Hs^{l- C)H + Ch, (5.4) 

in which C is an integration constant, and for simplification, and H2 are denoted as H 
and h. Using Eq. (5.4) we can from Eq. (5.2) obtain 

e = (1 - C)H'^ + 2Hh + Ch^ = KL, (5.5) 
K ^[{l-C)H + ahl (5.6) 
L = [H + bh], (5.7) 



where the constants a and h are the functions of C. Substituting the relation (5.5) into the 
viscosity function Eq. (2.11), then the Einstein's field equations (5.1) lead to 

'^^'^{l-C + a){KLYW -KW, (5.8) 

^ = -(l + 6)(i^L)^W- (5.9) 
dt 2 

Dividing Eq. (5.8) by Eq. (5.9), one gets 

dK a{l-C + a){KLY -2K 



dL a{l + h){KLY -2L 



(5.10) 
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After defining tlie variables 



A^L/B, (5.11a) 
S = (l + 6)X-(l-C + a)L, (5.116) 



Eq. (5.10) gives a simple form 

dA ail-C + a 



dB. (5.12) 



[^2 + (A/(l - C + a)]" 2^2"-^ 

The variables are now separated and the cqnation is Integratablc. (One can prove that 
1 — C + a is nonzero.) Using the methods described in the above sections we can therefore 
analyze any dimensional Bianchi type 1 cosmological models with energy density dependent 
bulk viscosity. The results show the same characteristics as those in the models with two 
Hubble functions. 



6 CONCLUSIONS 

We have analyzed in detail the anisotropic cosmological models with bulk viscosity {rj) which 
is a power-law dependence upon energy density (e), i.e., t] = a|e|", when the universe is filled 
with stiff mattery p = e. We are interested in the cosmological solutions that will eventually 
go to the states of positive Hubble functions in the latter stage. Although the exact solutions 
could be obtained only when the n is an integer, we are able to clarify the characteristics of 
evolution for the models of any n. 

Let us give a summary. (1) There have been two kinds of solutions in the n — model, 
which start either with diverge Hubble functions and infinite energy density at finite past 
or with finite Hubble functions and negative energy density at infinite past. However, both 
solutions are driven to a de Sitter space-time asymptotically. (2) All the solutions in the 
< n < 1/2 models will start with finite Hubble functions and negative energy density 
at infinite past and then are driven to a de Sitter space-time asymptotically. (3) For the 
1/2 < n < 1 models, the initial state is with finite Hubble functions and negative energy 
density at infinite past; however, they can go to the infinite expansion state, de Sitter space- 
time, or Priedmann universe at final stage. (4) The n > 1 models will always begin with 
Kasner-type curvature singularity at final past in which the energy density is zero, however; 
and then they are driven to the above-mentioned three kinds of states asymptotically. All the 
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solutions that begin with e j and then are attracted to a deSitter spacetime or Priedmann 
universe are free of cosmological singularity for any finite proper time. 

Historically, Murphy [8] presented the exact solution of the n = 1 model, and showed 
that the bulk viscosity can eliminate the big band singularity at any finite proper time. 
Belinskii and Khalatnikov [29] then analyzed the n = 1 Bianchi I model; they found that 
the cosmology is with the vanish energy density in the initial phase in which the Kasner 
singularity will arise. Now. the investigation in this paper shows that all the n > 1 Bianchi 
I models will share the same characteristics of initial singularity, and that the n < 1 Biancm 
models could give us some cosmological solutions that arc free of singularity for all finite 
proper time. However, these singular-free solutions have negative energy density in the early 
epoch. 

The models discussed in this paper arc only for the stitt matter As the case of stiff matter 
is special because the shear and matter density behave in the same way in the absence of 
viscosity and vacuum and nonvacuum perfect fluid solutions are formally similar, the same 
models, while with other matter fields are certainly interesting, remain to be studied. 

Finally, we want to mention that the prescription adopted in this paper can also be 
used to analyze the Bianchi type I cosmological models with energy density dependent shear 
viscosity [30] . The details will be presented elsewhere. 
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